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Abstract 

-t— > ■ 

d ■ We prove some new results which justify the use of interval truncation as a means of regu- 

larising a singular fourth order Sturm-Liouville problem near a singular endpoint. Of particular 
interest are the results in the so called lim-3 case, which has no analogue in second order singular 
problems. 

T- H ■ 

> 

^ ■ 1 Introduction 

■ In 1978, Bailey Gordon and Shampine |lj released a code (SLEIGN) for computing eigenvalues of 

Sturm-Liouville problems. This code was remarkable both for its reliability and for the fact that it 
was able to handle singular as well as regular problems, with a minimum of user input. SLEIGN's 
strategy for dealing with problems having singular endpoints was to truncate the interval near a 
singular endpoint, thereby regularising the problem. At the time, no rigorous proofs were given for 
the universal validity of this strategy. Recently, however, a new code (SLEIGN2) has appeared, 
capable of dealing with much more general singular second order problems: in particular, SLEIGN2 
can deal with non-Friedrichs boundary conditions near a lim-2 singular endpoint, and can deal 
with cases where there is an infinite sequence of eigenvalues tending to — oo. The development of 
SLEIGN2 followed the work of Bailey, Everitt, Weidmann and Zettl Q in proving rigorously the 
types of spectral convergence which could be expected from the interval truncation process near 
singular endpoints. 

In 1994, Greenberg and Marietta produced a code (SLEUTH) for solving fourth order regular 
Sturm-Liouville problems [ fl2|| . This was the culmination of three years of work on higher order 
self-adjoint ODEs Jl0|, |l) and numerical solution of eigenproblems for Hamiltonian systems H]. 
The new code gave sufficient increases in speed over the Hamiltonian systems code [Q for it to be 
feasible to solve singular fourth order problems; such problems were treated on an heuristic basis in 
[ p"2| . The purpose of the present work is to prove results similar to those produced by Bailey, Everitt, 
Weidmann and Zettl [0, but in the context of fourth order problems. For lim-2 and lim-4 singular 
endpoints we use methods which are direct adaptations to the fourth order case of the methods 
of However, for the lim-4 case, a new complication is the existence of complex as well as real 
boundary conditions. For the lim-3 case new difficulties are present, as this case does not arise for 
second order problems and cannot be treated by any of the methods described in [|| . We overcome 
these difficulties by using the oscillation theory described in [jllj . 
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tDept. of Mathematics, University of Leicester, University Road, Leicester LEI 7RH 
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2 ODE theory for the fourth order problem 

A fourth order Sturm-Liouville equation is an equation of the form 

£y = Xy, a<x<b, (2.1) 
where £ is a differential operator of the form 

*v ■= {[(p W)' - («(s)yT + • (2.2) 

Here 1/p, s, q and w are locally L 1 in (a, 6), with p and w positive almost everywhere, s and q 
real-valued. The endpoints a and b may be finite or infinite. The endpoint x = a is regular if it is 
finite and if there exists a' > a such that 1/p, s, q and w are in L 1 ^, a'); a similar definition holds 
for x = b. Any other sort of endpoint is called singular. 

Let L 2 (a, b; w) be the space of functions f(x) on (a, b) such that 

,6 

/ \f(x)\ 2 w(x) dx < oo. 

J a 

L 2 (a, b; w) is a Hilbert space with inner product 



(/,#) = / f(x)g(x)w(x)dx, 



and norm ||/|| = y/ (/, /). We shall say that a function / is square integrable at a if, for some 
P G (a, 6), 

|/(x)| 2 w(a;) < oo. 



Square integrability at 6 is defined similarly. 

Given a function y we define the quasi-derivatives 

= y, tf W = y M = w « y [3] = _(„,//)/ + sy '. (2 .3) 

These quasi-derivatives were introduced for scalar 2n th order problems by Naimark fl5| ; see also 
Everitt and Zettl @ and Zettl @ for further information on quasi-differential operators. 

Definition 2.1 The maximal domain D max is the set of functions y whose quasi-derivatives y^°\ 
y^\ j/ 2 l and y^ are all absolutely continuous, and for which y and ly lie in L 2 (a,b;w). The 
maximal operator L max is the operator defined by L max y = ty on the domain D(L max ) = D max . 

Definition 2.2 The pre-minimal domain C TO , n is the set of all functions in D max having compact 
support in (a, b) . The minimal domain D m i n is the closure of C m i n in the graph norm 

\\y\\G = (y,y) + (ey,£y)- 

The minimal operator L m i n is the restriction of L max to D m i n . 

It is not difficult to see that L m i n is the adjoint of L max . Every self-adjoint extension of L m j n is a 
restriction of L max to some domain D between D m i n and D max . It may happen that L max = L m i n , 
in which case both are self-adjoint and I has only this one self-adjoint realisation. However, L max 
is not generally self-adjoint. It is known that a self-adjoint extension L of L m i n has a domain 
D(L) C D max which is determined by certain boundary conditions at the endpoints. In order to 
know how many boundary conditions are required, we need the following classification of endpoints. 
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Definition 2.3 An endpoint is said to be of lim-p type if the space of solutions of the differential 
equation £y = iy which are square integrable at that endpoint has dimension p. 

For fourth order Sturm Liouville problems, the only possibilities are p = 2, p = 3, and p = 4. 
The lim-2 case is analogous to the limit-point case for second order equations; the lim-4 case is 
analogous to limit-circle. Lim-3 has no second order analogue. Regular endpoints are of lim-4 type. 
These endpoint types should not be confused with the deficiency indices of the minimal operator: 
for example, if both endpoints are of lim-2 type then the deficiency indices are zero. 

Boundary conditions are imposed by means of the Lagrangian form, which we now define. 

Definition 2.4 The Lagrangian form [/, g] of two functions f and g in D max is defined by 
[f,9](x) ■■= {/ [01 (^)F^ + / [11 W^)} - {f l3 \x)y^) + f^(x)' g m{xj}. 



with the notation of (2.5) for quasiderivatives 



Using another notation, 
where 



[f,g](x) = UfVg - VfUg, 



f[0] \ f f{3] 



and the superscript T indicates the transposed matrix. Integration by parts shows that for 
D max and a < a < (i < 6, 

r p — 

(£f)gwdx- f(e g )wdx=[f,g](f3)-[f,g}(a). (2.4) 

: J a. 

This shows that the limits 

[f,g](a):= ]im [f,g](x), [f,g}(b):= lim [f,g}(x) 

both exist. 

The domain D(L) of a self-adjoint extension of L m i n is determined by boundary conditions of 
the form [y, ip]{b) — [y,-0](a) = 0, where ip £ D max \D m i n . We are interested in separated boundary 
conditions, which are of the form [y, (f)](a) = or [y, ip](b) = 0. For such conditions, we may assume 
that (j) = in a left neighborhood of 6, and if) = in a right neighborhood of a. 

At a lim-2 endpoint no boundary conditions are required. At a lim-3 endpoint, one boundary 
condition is required. Suppose, for example, that x = a is a lim-3 endpoint. Then we can choose 
any (j) 6 D max \D m i n , such that [4>, 4>](a) = 0, and impose a condition 

M](a) = 0. (2-5) 

All valid boundary conditions at x = a would have this form. Now suppose that x — a were 
of lim-4 type. Then we would have to choose two maximal domain functions <j)\ and <f>2 with 
[<^i,<^i](a) = [4>i, <?y( a ) = [<t ) 2->4>2\{ a ) = 0, and such that no (non-trivial) linear combination of (f>\ 
and 4>2 is in D m i n . With these functions we would then impose two boundary conditions 

[y,ct> 1 }(a) = 0, [y,4> 2 ](a) = 0. (2.6) 



Equation (2.4) proves the following result which will be required later. 



Lemma 2.1 If yi and yi are two solutions of £y — Xy for the same real X then [t/i , j/2] is constant. 
If yi and y2 are solutions for the same complex X then [yi , j/2] "is constant. 
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We have now seen how to construct self- adjoint realisations L of I by restriction of L max to a 
domain D(L) C D max . As we mentioned earlier, self-adjoint realisations of i may also be obtained 
by extension of D m i n . It is particularly useful to view D{L) in this way when constructing a core 
for L. 

Definition 2.5 A set C C D(L) is said to be a core of L if D(L) is the closure of C in the graph 
norm. 

In essence there are only six possible combinations of endpoint singularities: lim-2-lim-2, lim-2- 
lim-3, lim-2-lim-4, lim-3-lim-3, lim-3-lim-4 and lim-4-lim-4. We shall consider each case in turn. 



Lemma 2.2 (Domains and cores of self-adjoint extensions) 

lim-2 lim-2 In this case D{L) = D min . Since D min is the closure of C„ 
C min is a core of L. 



in the graph norm, 



lim-2 lim-3 Suppose we have a boundary condition [y, "4>](b) = and suppose without loss of gen- 
erality that ip is zero on a right neighbourhood of x = a. Then D{L) = D min © Span(tp). The 
set C min © Span(ip) is a core of L. 

lim-2 lim-4 Suppose that we have boundary conditions [y, ipi](b) = = [y, ip2]{b) and suppose that 
ipi and ip2 we zero on a right neighbourhood of x — a. Then D(L) = D m i n © Span(ipi,ip2)- 
The set C m i n © Span(ipi,ip2) is a core of L. 

lim-3— lim-3 Suppose that we have boundary conditions [y, <j>](a) = and [y, ip](b) = and suppose 
that tp is zero on a right neighbourhood of x = a while <fi is zero on a left neighbourhood of 
x = b. Then D(L) = D min © Span(<f), ip). The set C min © Span((j), ip) is a core of L. 

lim-3— lim-4 Suppose that we have boundary conditions [y, <p](a) = and [y,tpi](b) = = [j/,-0 2 ](6) 
and suppose that <f) is zero on a left neighbourhood of x = b while tpi and ip2 are zero on a right 
neighbourhood of x = a. Then D{L) = D min ® Span(<f>, ipi , ^2) • The set C m i n © Span(<f), i/>i, 1P2) 
is a core of L. 

lim-4— lim-4 Suppose that we have boundary conditions [y, <j>i](a) = = [y, 4>2](a) and [y 7 ipi](b) = 
= [y, ip2](b) and suppose that <f>\ and (f>2 are zero on a left neighbourhood of x — b while ipi 
and ip2 are zero on a right neighbourhood of x = a. Then D(L) = D m i n © Span((f>i, (f>2, 4>i, ^2)- 
The set C min © Span(4>i, foipi, 1P2) is a core of L. 

We shall finish this section by describing some results from the theory of regular fourth order 
problems which we shall require later in this paper. We start with an alternative description of 
boundary conditions for a regular endpoint. Given a function y we define a vector z y by 



z y = (y[°],yW yflyPlf = { y iV > , - {py »)> + S y> iP y»)T . 

Recall that we have defined the vectors u v and v y by 

u y = (yM,yMf, v v = (yW,yMr. 

Define matrices J and 5* by 



(2.7) 



(2.8) 



J = 



( -1 



1 
\ 1 



-1 


J 







s = 



( Xw — q 


1 

\ 



0^ 

1 

l/p j 



(2.9) 
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The equation (2.1) is equivalent to the Hamiltonian system 

Jz' y = Sz y . (2.10) 

At a regular endpoint, say x — a, one imposes self-adjoint boundary conditions as follows. Let A\ 
and Ai be two 2x2 matrices such that the 2x4 matrix (A1A2) is of full rank (2) and such that 
A1A2 is Hermitian. Then the condition 

AtUyia) +A 2 v y (a) = (2.11) 

is a valid self-adjoint boundary condition at x = a, and every self-adjoint boundary condition at 
x = a has this form for suitable A\ and A2. 

Associated with the matrices A\ and A2 is a matrix solution Yl of the Hamiltonian system. Yl 
is a 4 x 2 matrix such that JY' L = SYl and 

Y L (a) = ("!!)■ (2-12) 
Usually Yl is partitioned into two 2x2 matrices Ul and Vl'- 



Y l = 



Ul 
Vl 



The matrix Ul is singular only at isolated points, and so we may define almost everywhere the 
matrix 

W L = YlUI 1 . 

It may be shown that Wl is Hermitian when A is real. Similar ideas hold for a right endpoint x = b, 
although it is then usual to denote the relevant matrices by Yr, Ur, Vr and Wr. 

In this paper we shall use extensively the oscillation theory developed in [jiof and summarised, for 
the case of fourth order problems, in pl[ . The theory was developed for real boundary conditions in 
[fiol , but it can be shown that the same proofs extend to the case of complex boundary conditions. 
The main result we shall require is the following. 

Theorem 2.1 Consider a regular fourth order problem over a finite interval (a, b) with separated 
self-adjoint boundary conditions at a and b. Let 



Sl(x,X)— {rank deficiency of U f{t , A)} . 

a<i<x 

Sr(x, A) = {rank deficiency o/[/^(<, A)} . 



x<t<b 

Then 



Sl and Sr are finite; 



• given any point c G [a, b] there is an integer a(c, A) whose value is either 0, 1 or 2, such that 
the number of eigenvalues of the eigenproblem strictly less than A is 

N(\) = S L {c, A) + 5 R (c, A) + a(c, A); (2.13) 

• if Ul(c, X) and Ur(c,X) are nonsingular then 

a{c, \) = v#(W l (c,\)-Wr(c,\)), 
where v#{W) = the number of negative eigenvalues ofW; 
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• if c = b and Dirichlet boundary conditions y(b) = = y'(b) are imposed then Sr = a = in 

• if c = a and Dirichlet boundary conditions y{a) = = y'{a) are imposed then 5l = cr = in 

fluty . 

Formulae for <r(c, A) covering all the other cases are given in |TT| . 

3 Operator convergence and spectral inclusion 

In this section we shall review some functional analytic concepts and results which will be required 
in later sections. This will be followed by some new results on spectral inclusion for singular fourth 
order operators approximated by regular fourth order operators. 

3.1 Functional analysis 

Notation We consider a self-adjoint operator L on a domain D(L) in a Hilbert space H with inner 
product (•,•). We also consider a sequence of self adjoint operators (Lj)°°^ 1 on domains D(Lj) in 
H. We shall denote by Sp(L) the spectrum of L and by Sp(Lj) the spectrum of Lj for each j. The 
notation || • ||g will denote the graph norm associated with L: 

\\y\\ 2 G :=(y,y) + (Ly,Ly), y e D(L). 



Definition 3.1 The sequence (Lj)j° =1 is strong resolvent convergent (SRC) to L if, for z G C\R 

and f G H , 

{L -zI)- x f -^(L-zI)-\f asj^w. 

A sufficient condition for strong resolvent convergence is given by the following result. 

Lemma 3.1 Suppose that C is a core of L such that, for each f G C, there is an integer N such 
that f G D(Lj) for all j > N , and 

lim Ljf = Lf. 

Then the sequence (Lj)°° =1 is SRC to L. 

See Reed and Simon @, Theorem VIII. 25] for a proof. 

Definition 3.2 We say that (Lj)jZ 1 is norm resolvent convergent (NRC) to L if, for some 
z G C\R, 

lim \\{Lj - zl)- 1 — (L — zl)- x \\ = 0. 

Definition 3.3 

(i) The sequence (Lj)j^_ 1 is spectrally inclusive for L if, for every A G Sp(L), there is a sequence 

with \U) g Sp{Lj) such that lim^oQ = A. 

(ii) The sequence (Lj)°° =1 is spectrally exact for L if it is spectrally inclusive for L and if, given 

a sequence (A^)^^ with X^' G Sp(Lj), every limit point of the sequence is in Sp(L). 

The reason for introducing the concept of SRC lies in the following result (see Reed and Simon 
fH Theorem VIII.24]). 
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Theorem 3.1 Suppose that L is a self-adjoint operator on a Hilbert space H and let Lj be a sequence 
of self-adjoint operators on H which is SRC to L. Then Lj is a spectrally inclusive sequence for L. 

Suppose also that we denote by -P(Ci)C2) the spectral projection of L associated with an interval 
(Ci ? C2) whose endpoints are not in the spectrum of L and let Pj(Ci,Cz) be the corresponding spectral 
projection for Lj, for each j. Then for any f G H, 

lim P J (Ci,Ca)/ = i , (Ci,Ca)/. 

3—>oo 

In particular, this implies that if X is an isolated eigenvalue of multiplicity k for L then there must 
be k eigenvalues (counted according to multiplicity) of Lj which converge to X as j — > 00. 

There is a corresponding connection between NRC and spectrally exact convergence, given by 
the following result (see Reed and Simon 0, Theorem VIII.23]). 

Theorem 3.2 Suppose that L is a self-adjoint operator on a Hilbert space H and let Lj be a sequence 
of self-adjoint operators on H which is NRC to L. Then Lj is a spectrally exact sequence for L. 



Furthermore, using the notation of Theorem 3.1 for the spectral projections, if (Ci 5 C2 ) * s an 



interval whose endpoints are not in the point spectrum of L then 

lim ||P J (Ci,C2)-P(Ci,C 2 )|| =0 



Remark . It may happen that L is unbounded below but that each Lj is bounded below. In this 
case it will be possible to index the eigenvalues of each Lj as, say, 



A§° < X\ J > - ■ 



assuming that these all exist. In this case, Theorem 3.2 would tell us that for each fixed k, 



lim 



U) 



-00. 



3.2 Fourth order Sturm-Liouville operators: interval truncation and spec- 
tral inclusion 

We construct regular approximations to our singular Sturm-Liouville problems by interval trunca- 
tion. We shall require two technical lemmas to allow us to set up approximating boundary conditions 
near lim-3 and lim-4 endpoints. 

Lemma 3.2 Suppose that the endpoint x = a is of lim-3 type for the differential expression £. Let 
4>{x) be a function in D max \ D m i n such that [cf>, 4>](a) = 0. Then there is a real function z{x) such 
that the boundary condition [y,(j)](a) = is equivalent to [y,z](a) = 0. A similar result holds for the 
endpoint x = b. 

Proof: We may suppose that <j>(x) = in a left- neighborhood of x = b. Let (j){x) = u(x) + iv(x), 
where u(x) and v(x) are real. Then 

[(f), <p](a) = [u, u](a) + [v, v](a) + 2i[v, u] (a) = 0. 

Therefore [v,u](a) — 0. Furthermore, 

[u,u](a) = = [v,v](a), 
since u and v are real. This implies that 

[u, 4>](a) = [u, u]{a) — i[u, v](a) = 0, 
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[v,4>](a) = [v,u](a) -i[v,v](a) = 0. 
Thus u and v satisfy the boundary condition [y, 4>](a) = 0, and therefore 

u = yi + ci<f>, v = y 2 + c 2 <j>, 

where 2/1,2/2 <= D min and Ci,c 2 are constants. If c\ = c 2 = 0, then both u and v E D min , and so 
4> E Dmin- Therefore one of the constants, say c\ 7^ 0. This implies that u E D max \ D m i n and the 
condition [y, <j>](a) — is equivalent to [y,u](a) = 0. □ 

Lemma 3.3 Suppose that the endpoint x — a is of lim-4 type for the differential expression I and 
suppose that boundary conditions \y,<f>i]{a) = = [y, <j> 2 ] {a) are imposed at x = a. Then for each 
real X there are solutions of the equation ty = Xy, say z\ and z 2 , such that the boundary conditions 
at x — a are equivalent to the conditions \y,Z\\(a) = = [y, z 2 ] (a). Similar results hold for the 
endpoint x = b. 

Proof: As we are interested in separated boundary conditions we do not want to be concerned with 
a possible singular endpoint at x = b; to this end we assume that x = b is regular, with two regular 
self-adjoint boundary conditions imposed there. If this is not the case we can always arrange for it 
to be so, by interval truncation, without making any changes to the boundary conditions at x = a. 
Consider the symmetric closed operator L a whose domain D a is the set of functions y E D max such 
that 

• y satisfies the regular boundary conditions at b; 

• [y,g]{a) = for an y a in A, 



■'max- 



Note that L a is a symmetric operator. For if 2/1,2/2 € D a , then 

rb rb 

/ {iyiW 2 wdx- I yi(£y 2 )wdx = [2/1,2/2] 0) = 0, 

since 2/1 and y 2 satisfy self-adjoint boundary conditions at x = b. The adjoint L* a of L a is the 
maximal operator associated with L a and its domain is the set of functions y E D max such that y 
satisfies the regular boundary conditions at b. 

The deficiency index of L a is the dimension of N(L* a — il). It is clearly two, because the fact that 
x = a is of lim-4 type means that the space of solutions which are square integrable at x = a is of 
dimension four; N(L* — il) is the two-dimensional subspace of these satisfying the two boundary 
conditions at x = b. We now fix any real A and any point c in (a, b) and construct four solutions 
2/1, . . . , 2/4 spanning the whole solution space for ly = Xy, say by using the initial conditions 

/ Vj(c) \ 

y'M 

-(py»y(c) + sy>(c) 
V *<(c) J 



where the ej are the standard unit vectors in R . For each i and j the Lagrangian form [yi, yj](. 
is constant, so it is easy to see that the matrix W — ([2/i, 2/j]) i s nonsingular. We now form four 
more functions 2/1, . . . ,2/4 which are functions in D max with yj(x) — yj(x) in a right neighbourhood 
of x = a (which includes the point x = c), and jjj(x) = in a left neighbourhood of x = b. The 
matrix of Lagrangian forms for these functions, W, has the property that 

W(b) - W(a) = -W(a) = -W(c), 
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and W(c) is nonsingular. By Lemma 10.2.17 of Hutson and Pym |L3[ this means that y\, . . . , 2/4 are 
linearly independent relative to D a , and by Theorem 10.2.18 of ]l3fl any two self-adjoint boundary 
conditions at x = a are equivalent to boundary conditions of the form 

[y,z 1 ](a)=0, [y,z 2 ](a)=0, (3.1) 

where z±, Z2 G Span(yi, . . . , 2/4). In particular, our boundary conditions [y, 4>\](a) — and [y, fa] (a) = 



can be replaced by conditions of the form (3T). Finally since ijj — yj near a for j = 1, . . . , 4, the 
proof is complete. □ 
Remark . The functions z\ , z-i in the previous lemma need not be real. There may be no real 
functions satisfying the lemma. 

We are now ready to construct regular truncated interval approximations to our singular prob- 
lems. We approximate the interval (a, b) by intervals (aj,bj), where a,j — > a+ and bj — > 6— as 
j — > 00. If x = a is regular then we can take a,j = a for all j, and if x = b is regular we can take 
bj = b for all b. Otherwise we must assume that Oj G (a, b) and bj £ (a, 6) with a-,- < bj for all j. 
At a regular endpoint the boundary conditions are inherited from the original problem. The list 
below describes what happens at singular endpoints. In each case, we take for a core of L the core 



described in Lemma 2.2 



lim-2 lim-2 In this case we impose boundary conditions y{aj) = = y'(aj) and y(bj) = = y'{bj). 
Any function in the core of L will have compact support in (a, b) and will therefore satisfy 
these boundary conditions for all sufficiently large j. 

lim-2 lim-3 Suppose we have a boundary condition [y,ip] (b) = where ip is zero on a right neigh- 



bourhood of x = a. By Lemma B.2 we may suppose that tfj is a real function, which implies 
that [ip, 4>](bj) = 0. Three of our boundary conditions are clear: we impose y{aj) = = y'{aj) 
and [y,ip](bj) = 0. For the fourth boundary condition we impose any self adjoint condition at 
x = bj. Any function in the core of L will certainly satisfy the conditions at x = Oj for all 
sufficiently large j; also, since any function y from the core is a multiple of near b it will 
satisfy [y, ip]{bj) = 0. Finally since the remaining boundary condition is a valid self-adjoint 
boundary condition, it is necessarily satisfied by ip, hence by multiples of ip, hence by y for all 
sufficiently large j. 

lim-2 lim-4 Suppose that we have boundary conditions [y, ipi](b) = = [y, fa]{b) and suppose that 
ipi and ip2 are zero on a right neighbourhood of x = a. The boundary conditions at x = cij 
are as in the first two cases. At x = bj we would like to impose the conditions [y, ipi](Jbj) = 
and [y,fa](bj) = 0. We know that [?/>i, fa](b) = but we do not know that [ip±, tp2](bj) = 0, 
so unfortunately the obvious boundary conditions may not be self-adjoint. To get round this 



we invoke Lemma 3.3, which tells us that we can assume that ipi an d ^1 are solutions of 



the differential equation with real A near 6, so that [ipi,ipi], [^1,^2] and [ip2 , 1P2] are constant 



(and therefore 0) near b by Lemma 2.1. Thus our boundary conditions [y, ^>i]{bj) — and 
[y, 1P2] (bj) — are self-adjoint. If a function y is in the core of L then near x = b it is a linear 
combination of ipi and ^2 and therefore satisfies the boundary conditions. 

lim-3— lim-3 Suppose that we have boundary conditions [y, (f>](a) — and [y, ip](b) — and suppose 
that ip is zero on a right neighbourhood of x = a while <f> is zero on a left neighbourhood 
of x = b. We may assume that <fi and ip are real. Then we impose boundary conditions 
[y, (j)](aj) — and [y,ip](bj) — 0, plus one other self-adjoint condition at each of Oj, bj. Any 
function in the core of L is a multiple of (j) near a and a multiple of ip near b and therefore 
satisfies the boundary conditions at aj and bj for all sufficiently large j. 

lim-3— lim-4 Suppose that we have boundary conditions [y, (p](a) = and [y, ipi](b) = = [y, 4>2]{b) 
and suppose that (f> is zero on a left neighbourhood of x — b while ipi and ^2 are zero on 
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a right neighbourhood of a; = a. We may assume that <j> is real. We impose the condition 
[y, 4>](a,j) = together with one other self-adjoint condition at aj. For the conditions at bj we 
assume once more that ipi and ip2 are solutions of the differential equation with real A near b, 
and we impose the conditions [y, ipi]{bj) = = [y,ip2][bj). Once more we can show that any 
function in the core satisfies the boundary conditions for all sufficiently large j. 

lim-4 lim-4 Suppose that we have boundary conditions [y,(j>i](a) = = [y,<j>2[(a) and [y,ipi](b) — 
= [y,ip2]{b) and suppose that <f>i and <p2 are zero on a left neighbourhood of x — b while ip\ 
and ip2 are zero on a right neighbourhood of x — a. We can assume also that <pi and </>2 are 
solutions of the differential equation with real A near x = a and that and ip2 are solutions 
near x = b, and then impose the conditions [y,<fri](a.j) = = [y,<j>2[{a,j) and [y,ipi](bj) = = 
[y, ij)2\(bj) at aj and bj respectively. Any function in the core of L will satisfy these conditions 
for all sufficiently large j. 

Of course, the problem of maintaining self-adjointness of the boundary conditions when truncat- 
ing at a lim-4 endpoint do not arise if the lim-4 endpoint happens to be regular, because then no 
truncation is required. 

Associated with each truncated interval problem is a self-adjoint operator Lj whose domain 
contains functions defined over [aj , bj] and satisfying appropriate boundary conditions at aj and bj . 
If we are to invoke the results of the previous section on SRC to obtain spectral inclusion, then we 
must create from the operators Lj new self-adjoint operators Lj which act on functions defined over 
the whole of (a, b). This is an uninspiring technical process which we shall now discuss. 

We shall consider L 2 (aj, bj; w) to be a closed subspace of L 2 (a, b; w) by extending the functions 
in L 2 (aj, bj; w) to be zero in (a, b)\(aj, bj). Then we have a splitting 



where 



L 2 (a, b; w) = L 2 (aj, bj; w) © L 2 (a,j, bj; w) , 



L (aj,bj;w) =L (a,aj;w)®L (bj,b;w). 
The operator Lj is extended to an operator L'j on L 2 (a, b; w) by defining it to be zero on L 2 (aj ,bj;w) : 

//, /,, O.,. (3.2) 

where Oj is the zero operator on L 2 (aj, bj; w) 1 - . The domain of L'j is 

D(L'j) =D(Lj)®L 2 (a j ,b j ;w) ± . 

It is easy to check that the L'j are self-adjoint. Unfortunately they also all possess as an eigenvalue 
of infinite multiplicity. This is no problem to us in the interpretation of the spectral inclusion results 
so long as we assume - as we may - that is not an eigenvalue of L. 

Lemma 3.4 In each of the cases described above, the operators L', are SRC to L. 



Proof: We shall use Lemma 3.1. Suppose that / lies in the core of L as described earlier. Our 
construction of the Lj has ensured that / satisfies the boundary conditions at aj and bj for all 
sufficiently large j, and so / lies in D(Lj). Thus we can compute Lf — Ljf, and from ( p.2|) we 
clearly have 



(Tf-T f\M - / for a; £ (dj,bj) 

[ J jJK } ~ \ {lf){x) {orxe(a,b)\(aj,bj). 



Thus 

I|l/-l;./|| 2 = r w(x)\e.f\ 2 (x)dx+ f w (x)\ef\ 2 {x)dx. 

J a J bj 

Since / lies in D max , if must be square integrable near a and b and so the integrals on the right 



hand side tend to zero as j — * oo. By Lemma 3.1 this completes the proof. □ 
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Corollary 3.1 The truncated operators Lj constructed above are spectrally inclusive for L. 



Proof: Theorem 3.1 implies that the L'j are spectrally inclusive for L, and therefore this is true for 



the Lj also. □ 

4 The Friedrichs extension for lim-2, lim-3 and lim-4 

We consider in this section a problem having one or two singular endpoints. We assume that the 
minimal operator is bounded below, so that it possesses a Friedrichs extension. Our objective is to 
approximate the eigenvalues below the essential spectrum of the Friedrichs extension. 

4.1 Friedrichs boundary conditions at every singular endpoint 

Let us start by setting up the appropriate domains for the minimal operators. Consider first the 
case of one singular endpoint, say x — b. We assume that two regular boundary conditions are given 
in the usual way at x = a. For our pre-minimal domain we take the set C m i n of all functions in the 
maximal domain which satisfy the boundary conditions at x — a and have compact support in [a, b). 
The minimal domain D m i n is the closure of C m in in the graph norm. For the case of two singular 
endpoints, we take the pre-minimal domain to be the set Cmin of all functions with compact support 
in (a, 6); the minimal domain D m i n is, once more, the closure of C m i n in the graph norm. In cither 
case we denote by L m i n the operator given on Drain by L m i n y = £y. Since L m i n is bounded below 
we may assume, by making a shift of q(x) in the expression for I if necessary, that 

(£y,y) > (y,y) \/y G D min . (4.1) 

We can then define the energy norm || • \\p by 

\\v\\E = (ty,y), y&D mm . (4.2) 

The energy norm is stronger than the weighted L 2 norm but not as strong as the graph norm. Recall 
that the Friedrichs extension is defined to be the operator whose domain is the closure of D min in 
the energy norm (see, e.g., Dunford and Schwartz 0, pp. 1240-1241]). We shall denote the Friedrichs 
extension by Lp and its domain by Dp . Obviously the energy norm defined in ( |4.2| ) extends to Dp. 

We now set up approximating regular problems. In the case of one singular endpoint, which we 
have taken for convenience to be x = &, we choose a monotone increasing sequence of points bj such 
that a < bj < b for all j and bj — ► b as j ' — ► oo. We set up regular problems on the truncated 
intervals [a, bj] which inherit their regular boundary conditions at x — a from the original problem 
and which have Dirichlet boundary conditions y = y' = at x = bj. For two singular endpoints we 
require a second (monotone decreasing) sequence of points aj such that a < Oj < bj < b for all j 

and aj — > a as j — > oo. In this case we also impose Dirichlet conditions y = y' = at the aj. We 

(i) 

denote by Lj the resulting regular operators and by their eigenvalues, for k = 0, 1,2, . . .; the 
eigenvalues of Lp will be denoted by when they exist. 

Finally, as a notational convenience, we shall define sets Hj as follows. In the case of one 
singular endpoint at x = b, Hj will be the set of maximal domain functions on [a, b) which satisfy 
the boundary conditions at x = a and have compact support in [a,bj); in the case of two singular 
endpoints Hj will be the set of maximal domain functions having compact support in (aj, bj). 

The following two lemmas are well known. The first is a standard result from the theory of 
regular Sturm-Liouville problems; the second comes from Berkowitz J^, Theorem 2.2] and the fact 
that Dp is the completion of C m i n under the energy norm. 

Lemma 4.1 For each k and j, 

(i) • c / W v \\e 

dim(v)=fe+i, vcHj Ivev \\ v \\ 
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Lemma 4.2 Suppose that Lp possesses a kth eigenvalue fXk strictly below any essential spectrum. 
Then 

fi k = mf <^ sup — =-= — — rj— 

dim(v)=fc+i, vcc min Uev \\ v \\ 

We can now prove the following useful result. 
Lemma 4.3 For each j and k, 



Moreover if fj, k exists and is below the essential spectrum of Lp then 

Mfe < K } 

for all j and k. 

Proof: Since bj is monotone increasing (and a,j, when required, is monotone decreasing) it is clear 
that Hj C Hj + i for all j, from which the first inequality holds from Lemma [4.1| . The second 
inequality comes from Lemma L2 and the observation that Hj C C m in for all j. □ 

We can now state and prove our main result. 

Theorem 4.1 Suppose that Lp possesses a kth eigenvalue fi k strictly below any essential spectrum. 
Then 

(i) 

lim fii = fj, k . 

Proof: From Lemma 4.2, given e > there exists a (k + l)-dimcnsional space V k C C m i n such that 

Hk > sup -ji-jj- - e. (4.3) 

Because V k e is finite-dimensional there exists a positive integer J k such that for all j > J|, V k e C -Hj . 
Thus from (fO), 

> mf <^ sup -— 

dim(y)=fe+i, vcHj Ivev \\v\ 

for all j > J k . The right hand side of this inequality is just fj,u — e, so 

Mfe > Mfe ' - e 

for all j > J k . Combining this with the second inequality in Lemma |4.3| yields the required result. 
□ 



4.2 Two singular endpoints with Priedrichs boundary conditions at one 
endpoint only 

Theorem |4.l| deals with two cases: the case of one singular endpoint, and the case of two singular 
endpoints where we want the extension corresponding to the Friedrichs boundary conditions at 
both ends. When we come to considering problems with one lim-2 endpoint and one other singular 
(either lim-3 or lim-4) endpoint, we shall find it useful to have a result on eigenvalue convergence 
in which we approximate a problem with two singular ends by a problem with one singular end. 
We shall therefore mention the case where a problem has one endpoint, say x = a, at which we 
impose Friedrichs boundary conditions (if required, i.e. if the endpoint is of lim-3 or lim-4 type), 
and where the remaining endpoint x = b is singular, but where no truncation is effected there. The 
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point x = b may be of lim-2 type, in which case no boundary conditions are required there, or it 
may be of lim-3 or lim-4 type. In the lim-3 case there will be a boundary condition at x = b of the 
form [y,i/)](b) = 0; in the lim-4 case there will be two boundary conditions at x = b, of the form 
[y, ipi](b) = = [y, ip2](b). We denote by Lp the self-adjoint operator arising from the equation, 
the Friedrichs boundary conditions at x = a, when required, and any other requisite boundary 
conditions at x = b. The domain C m i n of our preminimal operator we now define as follows: 

lim-2 at b: j/ G C m i n if and only if y is a maximal domain function and there exists a G (a, b) for 
which y has support in (a, b). 

lim-3 at 6: t/G C m in if and only if y is a maximal domain function, [y,ip](b) = 0, and there exists 
a G (a, b) for which y has support in (a, b). 

lim-4 at 6: y G C m i n if and only if y is a maximal domain function, [y, ipi](b) = = [y, ip2\{b), and 
there exists a G (a, b) for which y has support in (a,b). 

The minimal domain is, as usual, the closure of C m i n in the graph norm. 

We choose a sequence of points aj — * a, with aj G (a, b) for each j, and we set up truncated 
interval eigenproblems on the intervals (cij , b) by imposing the regular Dirichlet boundary conditions 

y(aj) = = y'(aj). 

We do not truncate at 6; the boundary conditions at b, if required, remain unchanged, so our 
approximating Sturm-Liouville problems are now also singular. We denote the associated operators 
by Lj . To recover the results of the lemmas above we need to re-define the sets Hj . 

lim-2 at b: In this case we take Hj to be the set of maximal domain functions with support in 
(a,j,b). 

lim-3 at b: In this case we take Hj to be the set of maximal domain functions with support in 
(aj, b) such that [y,tp](b) = 0. 

lim-4 at b: In this case we take Hj to be the set of maximal domain functions with support in 
(<Xj,b) such that [y,ipi](b) = = [y, ip 2 ](b). 

Denote by the eigenvalues of Lj and by fik the eigenvalues of Lp. We now have the following 
result. 

Theorem 4.2 Suppose that the minimal operator is bounded below and suppose that /ik lies strictly 
below any essential spectrum of Lp. Then 

lim^'= W . (4.4) 

Before proving this theorem we require a technical lemma. 

Lemma 4.4 Let C(-) denote the essential spectrum of an operator. Then C{Lj) C C(Lp) for all j . 

Proof: This is a simple application of the operator splitting technique of Akhiezer and Glazman 
[[n p. 520]. Let Dj be the set of all functions / G D{Lp) with quasiderivatives /^( a i) = f^K a j) — 
p 2 \aj) — /I 3 l(a_y) = 0; and let Mj be the restriction of Lp to Dj. Let Aj = Dj n L 2 (a, aj-,w), 
Bj = Dj fl L 2 (aj,b;w), and let Aj, Bj be the restrictions of Lp to Aj, Bj, respectively. Then 
Dj = Aj © Bj and Mj = Aj © Bj. The truncated operator Lj is a self-adjoint extension of Bj. Let 
Kj be any self-adjoint extension of Aj. Then Kj ®Lj and Lp are two self-adjoint extensions of Mj. 
By §, Theorem 1, §105], 

C(L F ) = C(Kj © Lj) = C(Kj)UC(Lj). 
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Therefore C(Lj) CC(L F ). 



□ 



Proof of Theorem fi-'2\ Because /i^ lies below the essential spectrum of Lp, Lemma 4.2 still holds. 
Thus, given e > there exists a (k + l)-dimcnsional space V£ in the preminimal domain C m i n such 
that 

»k > sup Mf - e. (4.5) 
vev* \\v\\ 

Because V fe e is finite-dimensional and because each element of C m i n belongs to all the Hj for all 
sufficiently large j, there exists a positive integer J| such that for all j > J|, V£ C Hj. Thus from 



dim(v)=fc+i, vcHj 



sup ■ 



V\\E 



for all j > Jf.. Rearranging this result gives 



inf 



sup ■ 



\He 



(4.6) 



dim(V)=fc+i, VdH 3 {vev \\v 

We may assume that e is so small that \ik + e lies strictly below the essential spectrum of Lp and 
hence, by Lemma 4.4, below the essential spectra of all the Lj] in particular the left hand side of ( |4.6|) 
lies strictly below the essential spectrum of Lj. Standard variational theory (see, e.g. BerkowitzB 
Theorem 2.2]) now tells us that the eigenvalue exists and is given by the expression on the left 
hand side of (4.6), whence we deduce that 



for all j > Jf, . The inclusion Hj C C„ 
yield the inequality 

for all j > Jf:. This completes the proof. 



and the now established variational expression for ji^ 1 now 



□ 



5 The lim-4 case 

5.1 Truncation at the lim-4 endpoint 

We consider here the problem where the endpoint x = b is lim-4, and the other endpoint is lim-2 
or lim-3. We shall truncate only at the lim-4 endpoint x = b. Later we shall consider double 
truncations. L will denote a given self-adjoint extension of L m i n . 

Lemma 5.1 Let x = a be a lim-2 or lim-3 endpoint, and x — b a lim-4 endpoint. Let A G C \ R, 
and let L be a self-adjoint extension of L min . There are two independent solutions <f>\ (x) , fa (x) of 
the equation £y = Xy, such that: 

(1) fa, fa £ L 2 (a,b;w); 

(2) fa and fa satisfy the boundary conditions for L at x — a (if any); 

(3) [fa,fa](a) = 0. 

Proof: If the endpoint x = a is lim-2, there exist two independent solutions fa, fa € L 2 (a,b';w) 
for any b' G (a, b). These solutions also lie in L 2 (b' , b; w), since x = b is a lim-4 endpoint. Thus they 
lie in L (a, b;w). In this case there are no boundary conditions at x — a, and [0i,<?>2]( a ) = 0, since 
[yiiV2\ (a) = for all functions y x ,y 2 G D max . 
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Suppose that x = a is a lim-3 endpoint. By Lemma 3.2 the boundary condition at x = a is of the 
form [y, z](&) = 0, where z(x) is a real function. This implies that if 4>{x) satisfies the boundary 
condition, then <fi(x) does also. 

There are three independent solutions in L 2 (a,b;w). Two linearly independent solutions <fii,4>2 
satisfy the boundary condition at x = a. Since (f> 2 also satisfies the boundary condition, [<f>i, </> 2 ](a) = 
0. " " □ 



Lemma 5.2 Let x — b be a lim-4 endpoint, and let A 6 C. Suppose that <fi is a solution of iy = Xy 
such that [y, 4>]{b) = for all solutions y. Then tf> = 0. 

Proof: For any two solutions y%,y2, [yi, 2/2] ( x ) i s constant. This defines a skew-symmetric bilinear 
form on the space S of solutions. This form is nondegenerate in the sense that if z G S and [y, z] = 
for all y G S, then z = Q. The nondegeneracy can be checked as follows. Let a < c < b, and for 
j = 1, 2, 3,4, let yj be the solution of ly = Xy with the following quasi-derivatives (using the notation 
ofQ): 

(yf ] (c) yf(c) yf{c) yf{c)) T = e 3 , 
where {ei, 62, e3, 64} is the standard basis for C . The matrix A = ([yj, y k ]) is 



A = 



h 

-h 



where I2 is the 2x2 identity matrix. Since det A / 0, the bilinear form [y, z] is nondegenerate. Since 
[y,~4>] = for all y G S, it follows that = 0. □ 



Definition 5.1 Let x = b be a lim-4 endpoint, and let L be a self-adjoint extension of L m i n . 
The boundary conditions for L at x = b will be called real if they are equivalent to conditions 
[y, Zi](b) — = [y,Z2](b), where z\ and z 2 are real functions. Otherwise the boundary conditions will 
be called complex. 



Lemma 5.3 Let x — b be a lim-4 endpoint. Let X G C\R, and let L be a self-adjoint extension of L m y 

(1) There are two independent solutions Tp\(x), ip2(x) of the equation £y = Xy, such that 

(a) ipi, ip2 G L (a',b;w) for any a' G (a, b); 

(b) -01 and ip2 satisfy the boundary conditions for L at x = b. 

(2) [ipi, ip2\(b) — if and only if the boundary conditions at x — b are real. 

(3) If the boundary conditions at x = b are complex, thenipi, "02 can be chosen so that [^i,ip 2 ]{b) = 1. 

Proof: (1) Let S be the 4-dimensional space of solutions of £y — Xy. Since x — b is lim-4, all 
functions in S belong to L 2 (a', b; w) for any a' G (a, b). There is a 2-dimensional subspace T C S of 
functions satisfying the boundary conditions [y, z\](b) = = [y, z 2 ](6) of L at x = b. Take tpi, ip% to 
be any two independent functions in T. 

(2) If the boundary conditions at x — b are real, then ip satisfies these conditions if ijj does. Thus 
02 satisfies the boundary conditions, and so [-0i,0 2 ](^) = 0- 

Conversely, suppose that [ip 1 ,ip 2 ](b) = 0. Since [i/>2, 2 ](^) = e^so, 0> 2 satisfies the boundary 
conditions at x = b. Similarly, ip 1 satisfies these conditions. 
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We claim that the boundary conditions [y, Z\\{b) = = [y, Z2\{b) for L at x — b are equivalent to the 
conditions [y, (b) = = [y,ip 2 ](b). To show this, let a < a' < b. For a' < x < b, 



ipi(x) 



yi(x) + cizi(x) + c 2 z 2 (x), 
y 2 (x) + dxZi(x) + d 2 z 2 (x), 



(5.1) 
(5.2) 



where y\, y 2 S D m i n . By Lemma 5.2 a nontrivial linear combination of ipi and ip 2 cannot lie in D m i n . 
Therefore equations ( p.l| ) and (:k2) can be solved for z\ and z 2 in terms of "01, ip2 and functions in 
D m i n . This implies that the boundary conditions are equivalent to [y, ipi](b) = = [y,%^ 2 ]{b). 

Let ipi = u\ + ivi, ip2 — u 2 + iv 2 , where the functions Uk and Vk are real. Since i/j^, ip 2 satisfy the 
boundary conditions, the real functions u\, V\, u 2 , v 2 also satisfy the conditions. Two of these func- 
tions must be independent relative to D m i n . For otherwise ip-y and ip 2 would be dependent relative 
to D m i n ■ Therefore xj)\ and ip 2 can be expressed as linear combinations of two of these real functions 
modulo D m i n , and so the boundary conditions can be given by real functions. This proves (2). 

(3) If the boundary conditions at x = b are complex, then [ipi,ip 2 ](b) = a ^ 0. The function 
-01 = -j^>i satisfies [ijji,4> 2 ](b) = 1. Thus if ipi is replaced by ipi, then the required equation is 
satisfied. □ 



We will now calculate the Green's function for L. We shall first calculate it for real boundary 
conditions at x = b, and then indicate the result for complex boundary conditions. The equation 
ly — Xy = / has the form: 

[{py"Y - i s y')}' + (q- ^ w )y = fw. (5.3) 

This is transformed to Hamiltonian form as follows. Corresponding to the function y(x), consider 
the quasi-derivatives 

ui =y, «2 = v', v i = s y' - (py"Y, v 2 = py", 



Equation (5.3) is equivalent to 



Ul 
"2 



V2 



Jz = Sz + fw, 



where z = z y , J and S are as in (2.9), and 



/=(/,0,0,0) T . 



(5.4) 



Let 4>\, fa , ipi, ip2 be the solutions of ly = Xy from Lemmas 5.1 and 5.3. We obtain a fundamental 
matrix 

' U L U r 



$ = 



V L Vr 



where 



Ul = (u ( j }1 W0 2 ) 



31 92 
*'l <t>2 



Ur - {u^ lty a ) = I ^, ^ , . 



Vr = (v i>1 v^ 2 ) 



- (M7 - (p^y 



In the following, A T denotes the transpose of the matrix A 



- (p^'{y s v 2 - (p^'iy 
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Lemma 5.4 Suppose that the boundary conditions at x = b are real. Then: 

(1) ulv L - v£u L = o = ulv R - v£u R . 

(2) The functions ipi,ip2 can be chosen so that UrVl — VrUl — I (the identity matrix); 

(3) I/V'1,'02 are chosen to satisfy (2), then = I r] ff 



Proof: (1) By Lemma 2T, if and ip are solutions of £y = Ay, then [0, ip] (x) is constant. 

Since [0i,0 2 ](a) — = [i/>i , ^2] (&) 1 ^ follows that [0i, 2 ](£) = = [V>i, V^K^) f° r au ^ G ( a >^)- 

Furthermore, [0, 0](x) = for all functions 0(x). This implies (1). 

(2) Let S denote the 4-dimensional space of solutions of the equation ly — Ay. For £, rj G 5, 
[C,??](a;) = constant. We shall denote this constant by [C17]. The map (£,7 7) >— > [C,?7] defines a 
skew-symmetric, bilinear form on S x S. It was shown in the proof of Lemma 5.2 that this form is 
nondegenerate. 

Let S l = Span(0i,02), and Sr — Span (/0i, ip2 )■ If C £ £l H Sr, then £ is an eigenfunction 
with eigenvalue A for the self-adjoint operator L. Since A G C\R, this implies that £ = 0, and 

s L n s fl = {0}. 

If Ci V £ Sl, then [C,?j] = 0. Furthermore, is maximal with respect to this property: if £ G 5 
and [Cry] = for all t) £ Si, then ( 6 Sj. (This follows from the fact that if V is a maximal 
isotropic subspace of a symplectic space W, where dimVF = 2n, then dimV = n.) Sr is also 
maximal with respect to this property. This implies that [C,fj] defines a nondegenerate bilinear form 
on Sr. x Sl- Therefore, corresponding to the basis {01,02} for Sl, there is a dual basis {V'lW'f} 
for S R : [^,0,.] = % = 1,2). If we replace {^1,^2} by {V»f,Va}. then PflVt " Va^i = 7 - 

(3) From (1) and (2) it follows that 

~Vr UT\(U l Ur \ ( h 2 



vl -ul )\v L v R ) - \ o 2 j 2 ) - u - 

where /„ denotes the n x n identity matrix, and 0„ denotes the n x n zero matrix. □ 



We will now solve equation (5.4) by variation of parameters. Let 

ax) 

ly-(x) 



z(x) = *(x)[ it\ ), (5.5) 



—1 tuy ^<iu 



where 



are unknown functions. In the notation preceding equation ( p.4[ ), z — z y , where 

y = Ci0i + C202 + m^i + V2ip2- (5.6) 

Since y{x) must satisfy the boundary conditions at the endpoints, we require 

C(6) = = via). (5.7) 



Substituting ( |5.5| ) into (|5.4|), and using the fact that the columns of $ satisfy the corresponding 
homogeneous equation, we obtain 
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or 



rf{x) 



( °\ 



-/ 

V 0/ 



By Lemma 5.4, if the boundary conditions are real then 



T)'(x) 



( 






( -A 


\ 









-ifa 




-/ 






fa 






I 




{ <h 


) 



fw. 



Therefore, using (5.7), we have 



f(t)w(t)dt, r,(x) = 



Mt) 
Mt) 



f(t)w(t)dt. 



^From (5.6) we see that the solution y{x) of equation ( p.4| ) is 



G(x,t)f(t)w(t)dt, 



where the Green's function G(x, t) for real boundary conditions at x = b, is 

G(x,t) - 



4>i(x)ipi(t) + (f>2(x)ip2{t) for a < x < t < b 
<j>i(t)ipi(x) + 4> 2 {t)ip2{x) for a < t < x < b. 



(5.8) 



Letting S = (L — A) 1 , we see that for / 6 L 2 (a, 6; w), 



(Sf)(x) = (L — A)-V(.x) = / G(x, t)f(t)w(t) dt 



The calculation for complex boundary conditions at x — b is similar. Lemma 5.4 remains the same, 
except that 

UrV r - V^U R = K 



1 

-1 



and 

This implies that 

C(x) = 



Mt)-Mt) 



VI - KVl 

vl 



f(t)w(t) dt, n(x) = 



f(t)w(t) dt, 



and the Green's function for complex boundary conditions at x = b is 

G(x,t) -- 



Mx)^\{t) + 02 0*0^2 (*) + M x )Mt) - Mx)Mt) for a < x < t < b 
Mt)M x ) + Mt)M x ) for a < t < x < b. 



(5.9) 



We shall now truncate the interval near x = b. Suppose that the boundary conditions for L at x = b 
are [y, 0i](b) = = [y, 02](b), where 6\ and 9 2 are solutions of ly — Xoy for some real Ao- As before, 
we shall first carry out the calculation for real boundary conditions, and then indicate the result for 
complex boundary conditions. 

Let bj /* b, and consider the truncated operator Lj defined on (a, bj] with the same boundary 







[y,e 2 ]{bj) at x 



conditions as L at x = a, and the boundary conditions [y, 6i\(bj t 
The Green's function Gj(x,t) will resemble G(x,t) in (5.8), except that ipi and ip2 are replaced by 
solutions and if)^ which satisfy the boundary conditions for Lj at x = bj. 
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Lemma 5.5 Suppose that A G C\R. Then there are two linearly independent solutions of £y = Xy 
which satisfy the boundary conditions of Lj at x = bj and are of the form 



where 



= flh + <t>l + 4 ?) ^2 , ^ = 4>2 + 4 J) fa + fa , 



lim c[ 3) = lim cSp = lim d{ 3) = lim d^ = 0. 



(5.10) 



Proof: We shall solve for the constants cj!\ d ( k 3) in ( |5~To| ). Setting , 6i](bj) = = [V'F, 2 ](bj) 
and [V>2 , = = [V , 2 J \ ^2](&j), we obtain the following equations: 



i>i,0 2 ](M 



,0') 



= -fofc,02](&i), 
= -[^2^2](&i). 



There is no nontrivial solution to the equations 



H 



(5.11) 



c 1 [fa,9 1 ](b ] )+c 2 [fa,9 1 }(b 3 ) = 0, 
ci[<f>u6 2 \(b j ) + c 2 [cf) 2 ,6 2 \(b j ) = 0, 



since this would imply that c\4>\ + c 2 fa is an eigenfunction of the truncated problem with eigenvalue 
A e C \ R. Therefore 



A(&i) 



and similarly 



A(b) 



[<h,9i](b) [<t>MQ>) 
[fa,9 2 ] (b) [fa,o 2 ](b) 



^0. 



Therefore the equations (5.11) have a unique solution: 



Si) 



A(6i) 



.0) _ 



A(6j 



[fc,0i](6 3 -) 

[<feA] MM (bj) 

[V>iA](&i) [fa,6 x ](bj) 



and similar formulas for of and d^ ■ Since linx,-—^ [fa , 9k] (bj ) = [fa , 9k] (b) , and lim 
[ipi,6k](b) = 0, it follows that lim cf ] = lim c 2 j) 



J— >oo 



lim 4 j) 



lim d ( 2 3) = 0. 



bl>i,0k]Q>j) 



□ 



We now can find the Green's function for the truncated operator Lj . It follows the same pattern 
as (Bjh ; the formula is 



Gj{x,t) = G(x,t) + c\ 3 > 0i(a:)0i(t) + c< 3) </>i (x)0 2 (t) + 4" <fe 0)01 (t) + d% '> <t>i{x)4>2{t) for a < x < t < b 3 . 



G s {x,t) = G(x,t) + c< 3) 0i(t)0iO) + 4 3) 0i(t)02(z) + d[ j) Mt)M*) + 4 3 ' M*)Mx) 



for a<i<a;<6j. 



The above formula was calculated for the case of real boundary conditions at x = b. An analogous 
calculation shows that the formula for complex boundary conditions is the same, with constants 
that approach as j — * oo. 
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As described in section 3.2, we consider L 2 {a,bj\w) C L 2 (a,b;w), and the splitting 

L 2 (a,b;w) = L 2 {a, bj-,w) © L 2 (bj,b;w). 

The operator Lj with domain D{Lj) C L 2 (a,bj]w) is extended to L'j with domain D(Lj) 
L 2 (bj, b; w) by the formula 



3 

'r tt' \ ] w ■ ' I 1 i \ I / . "I ■ 1 1 i 1 1 1 - i 



where Oj is the zero operator on L 2 (bj,b;w). If P, denotes the projection of L 2 (a,b;w) onto 
L 2 (a,bj;w), then = Lj.Bj. Setting Sj = (Lj — A)^ 1 and Sj = (L'j — A) -1 , we have 



W - A )/ = { 

^(/) = ( j l;--a)- 1 / = { 



(Lj - A)/ for / e L 2 (a, w) 
-A/ for / e L 2 (bj,b; w), 



(L.-A)- 1 / for/eL 2 ( a ,6, ;U ;) 
-if for/eL 2 (6„6;u;), 



Pj<Sj — SjPj — / '; • 

Let A = (a, b) x (a, 6), A, = (a, 6j) x (a, &,) and set 

, ( + -,_{ Gj(x, t) for (x, t) e Aj 
^j{x,i) - | Q for ^ 6i ^. 

Then 

(S j P j f)(x)= [ G' 3 {x,t)f(t)w{t)dt. 

J a 

If L has real boundary conditions at x — b, then 

!c[ j) <t> 1 (x)</> 1 (t) + 4 J) 0l(^)02(t) + d[ ]) <t> 2 {x)<p 1 (t) + d l 2 j) <f> 2 (x)4> 2 (t) for a < x < t < bj 

c{ j) l p 1 (t)<j} 1 (x) + c ( 2 j) <j} 1 (t) l p 2 (x) + d[ 3) 4> 2 (t)4> 1 (x) + d l 2 i] 4> 2 (t)4> 2 (x) for a < t < x < bj 

<Ma;)iMt) + <l>2{x)i> 2 {t) ' for x < t, and bj < t < b 

^i(x)<t>i(t) + i> 2 (x)<fi 2 (t) for t < X, and bj < x < fc. 

The formula shows that G — G^e L 2 (A; w(x)w(t)) and G - G^- -> in the L 2 norm as j — ► oo. An 
analogous calculation shows that the same is true for complex boundary conditions at x = b. This 
proves the following theorem. 

Theorem 5.1 

(1) S — SjPj is a Hilbert- Schmidt integral operator. 

(2) S — SjPj — > in the Hilbert- Schmidt norm. 
The theorem implies the following corollary. 
Corollary 5.1 

(1) The sequence of truncated operators Lj is spectrally exact for L. 

(2) Lj and L have the same essential spectrum. 
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5.2 The lim-4, lim-4 case 

We now suppose that both endpoints are lim-4, and we shall consider double truncations. In this 
case, L has discrete spectrum. By translating the operator if necessary, we may suppose that is 
not an eigenvalue. The eigenvalues of L may be indexed by positive and negative indices: 

• • • < A_3 < A_2 < A_i < < Ai < A 2 < A 3 < • • • < \ n < ■ ■ ■ 

It is possible that L is not bounded below, and there are infinitely many negative eigenvalues. 
Let G(x, t) be the Green's function for L at A = 0. Thus 



(L-\f)(x) = / G(x,t)f(t)w(t)dt. 

J a 



There exist functions <p\, 4>2, ipi, tp2 G L 2 (a,b;w) which are solutions of ly = 0, and such that 
0i, <f>2 provide the boundary conditions for L at x = a, and ipi, ip2 provide the boundary conditions 
at x = b. G(x, t) can be expressed in terms of these functions by formulas similar to those in the 
previous section. There are now four cases, depending on whether the boundary conditions are real 
or complex at x = a and x = b. We shall not require the precise formulas. 

Let dj \ a, bj / b and consider the truncated operators Lj on the intervals [cij,bj] with 

boundary conditions [y, 0i](a.,-) = = [y, <p2]{aj) and [y, ipi](bj) = = [y, ip2](bj)- The solutions 

(i) (i) (i) 
<j>Y , <i>2 of £y — satisfying the boundary conditions at x = a,j are simply the restrictions 0- = 

4>i\[aj,bj]; and similarly the solutions satisfying the boundary conditions at x = bj are ip^ — 

Let Sl = span(0!,0 2 ) and Sji = span (ip\, ^2)- S'l H Sr = {0}, since is not an eigenvalue of 
L. This implies that is not an eigenvalue of Lj. Let A = (a, b) x (a, b) and Aj = (a,j, bj) x (a,j, bj). 
The Green's function Gj(x,t) for Lj at A = is defined in Aj and coincides with G(x,t) there. 

We have the splitting 



where 



Let Sj = Lj 1 and 



L 2 (a,b;w) = L 2 (aj,bj]w) L 2 (a J ,b j ;w) ± , 
L 2 (aj 1 bj]w) ± = L 2 (a, a,j\w) © L 2 (bj,b; w). 



Sj — Sj (B 



where Qj is the zero operator on L 2 (a,j,bj:w) ± . Let Pj be the projection of L 2 (a,b;w) onto 
L 2 (aj,bj-,w). Then PjS'j = SjPj = SjPj. SjPj is an integral operator with kernel 



G' 3 (x,t) = I G( p'^ f ° r G A -" 



for (x,t) G A\ Aj. 

It is now clear that G — G'j G L 2 (A:w(x)w(t)) and G'j — > G in the L 2 norm. We have proved the 
following. 

Theorem 5.2 SjPj — > S in the Hilbert- Schmidt norm. 
Let the eigenvalues of Lj be: 

• • • < A°] < < \ u \ < < A^ < \ { 2 j) < X ( 3 j) < ■ ■ ■< A« < • • • 

Since Lj is a regular operator, it has only finitely many negative eigenvalues. Nevertheless, for 
any (positive) integer fc, if A_fc exists, then for sufficiently large j, A^| exists, and X^P k — > A_fe as 
j 00. 
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(i) 

Theorem 5.3 For k > 0, A^ ; -> A fe as j — > cxd. Furthermore, if X-k exists, then for sufficiently 
large j, Lj has an eigenvalue A^L and A_jj. — > A_fe as j — > oo. 



Proof: 5 and are Hilbert-Schmidt operators, and SjPj — > 5 in the Hilbert-Schmidt norm. 

These are self-adjoint, compact operators with eigenvalues ^ = 1/Afc and /z^ = 1/Ajj. , respec- 
tively. (SjPj also has as an eigenvalue of infinite multiplicity, but this is not related to Lj or L, 
and can be ignored.) The negative eigenvalues are given by a min-max variational principle, and 
the positive eigenvalues by a max-min principle. For k > 0: 

/ifc = max min (Sx, x), 

Ik ll=i 

/i_fc = min max(Sa;, x), 

M k 

Ik ll = i 

where runs through the fc-dimensional subspaces of L 2 (a,b;w). The ^ are given similarly in 
terms of (SjPjX,x). The result now follows from the fact that 

\((S-S 3 P 3 )x,x)\ < 115-5^-11-0, asj-oo, 

for Ibll = 1. □ 



6 One singular endpoint: the lim-2 case 

In this section we consider a problem with one regular endpoint and one lim-2 singular endpoint. 
At the regular endpoint, say x = a, we impose two regular self-adjoint boundary conditions. The 
pre- minimal domain is then precisely the set C m i n which we described in Section ^ for the case of one 
singular endpoint, and the minimal domain is its graph-norm closure. Because the endpoint x = b 
is lim-2 the minimal operator possesses only one self-adjoint extension L. If the minimal operator 
is bounded below, then L is the Friedrichs extension, and the results of Section |] show how we 
can obtain spectral exactness below the essential spectrum; otherwise there is spectrum extending 
to — oo, possibly with gaps. If there are no gaps then spectral inclusion evidently implies spectral 
exactness; in the event of gaps, we will show that spectral exactness may not be obtained. 

The regular approximations Lj to L can still be formed as in Section ^. Any element of the 
pre-minimal domain - which is a core for L in the lim-2 case - will satisfy the Dirichlet conditions 
yftj) = y'ibj) = for all sufficiently large j, and so its restriction to the intervals [a,bj] will be in 
the domain of Lj for all sufficiently large j. This means that the spectra of the Lj give spectrally 
inclusive approximations to the spectrum of L. In particular, we have the following result. 

Proposition 6.1 Suppose that the minimal operator is unbounded below and, for each fixed integer 
k>0, let /ijL be the kth eigenvalue of Lj. Then 

lim iir' = — oo. 

It is also not difficult to see that if there are gaps in the essential spectrum of L then the bj may be 
chosen to ensure that for some point A* in one of these gaps, there is some k such that 



To see how this may be achieved let /3 be less than b and let fik(P) be the fcth eigenvalue for the 
problem on [a, j3\ with Dirichlet conditions at [3. By taking k sufficiently large we can ensure that 
/Ufc(/3) > A*; with k now fixed we can let (3 increase towards b until /ifc(/3) attains the value A*. Thus 
we obtain the following result. 
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Proposition 6.2 Suppose that L is unbounded below and suppose that A* is not a spectral point of 
L. Then the approximating operators Lj may be constructed so that A* is a spectral point of every 
Lj. 

The approximation of essential spectrum evidently requires something more than eigenvalues: 
it requires the spectral function or the Titchmarsh-Weyl M(A) matrix. Computing the spectral 
function is a difficult problem even in the second order case. For more information about ways of 
computing the M(A) matrix in the fourth order case, see Q. 



7 One singular endpoint: the lim-3 case 

We shall now consider a problem where the endpoint x = o is regular, and x = b is lim-3 singular. 
A self-adjoint extension L of L m i n will have boundary conditions of the form 

Aiu y (a) + A 2 Vy(a) = 0, (7.1) 

[y^](b) = 0, (7.2) 
where ip is a real function in D max \ D m i n . In the following, we will need to compare L with another 



self-adjoint extension Lp. Let Lq be the extension of L m i n with boundary conditions (7.1) at x = a. 
If L is bounded below, then so is Lq, and therefore Lq has a Friedrichs extension Lp. This operator 
will play an important role in the following. 

Suppose that L is bounded below, with a number of eigenvalues strictly below the essential 
spectrum - say Ao, . . . , A„. Note: we are not assuming that these are the only eigenvalues below 
the essential spectrum. We shall construct a sequence of regular operators Lj on truncated intervals 
[a, bj), such that for each k between and n, the fcth eigenvalue of Lj converges to the fcth eigenvalue 
of L as bj y b. This will require a careful choice of boundary conditions at the endpoints bj, which 
we now describe. 

We shall denote by 

' U L (x) ' 



Vl(x) 

a 4 x 2 fundamental matrix solution of the Hamiltonian form of the Sturm-Liouville equation which 
satisfies the initial conditions ( 2.1 2| ) at x = a. The dependence on A is suppressed in this notation. 



Associated with this fundamental matrix is the matrix 

W L (x) = V L {x)UZ\x), (7.3) 

which is defined except at a finite number of points in (a, b) (assuming that A lies below the essential 
spectrum). Suppose that, in the Hamiltonian formulation, the boundary conditions for Lj at x = bj 
have the form 

v = W R u, (7.4) 

where Wr is a real, symmetric matrix: 

W R =[l J). (7.5) 

By Theorem 2.1, the number of eigenvalues of Lj which are strictly less than A is 

N(b j ,X) = 5 L (b j ,X)+a(b j ,X), (7.6) 

where 

SL.(bj, A) = (Rank deficiency of Ul{x))- 

a<x<bj 
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If det U L (bj,X) # 0, then 

<r(b j ,\) = u#(WL(b j )-W R ) ) 

where v#(W) is the number of negative eigenvalues of W. We shall investigate the possible values for 
cr(bj, A) corresponding to different choices of the matrix Wr which defines the boundary condition 
at x = bj. In particular, we shall show that a — 1 is always possible. We shall use the notation of 
( [2.3] ) for quasiderivatives. We require the vectors 

V ,[0] \ / ^[3] 

V ,W J > v i' - I ^[2] 



Because the boundary condition defined by (7.4) must include the boundary condition defined by 
(ff^), we must have 

«ty = W R u^. (7.7) 

If we assume that ?/>M is non-zero, then we can rearrange this formula to obtain \x and v in terms of 
the unknown k and the known V , ' 1 ' J 4^ an d V^ 3 ' : 

„ = (^ - K^)/^\ V = |^ - ^(^ - ^ [01 ). 

To determine the number of negative eigenvalues of Wl — Wr we require the trace and determinant 
of this matrix. Expressing these in terms of n we obtain 

traced - Wfc) = trace(^) + - k (l + g^J . (7.8) 

det(W £ - W R ) = C + — ^-(^ - u^Witty). (7.9) 
The constant C in ( |7.9| ) is given by 

C = dct{W L ) + (fc(^ [0 V [3] - V [ V 2] ) - W [31 ) 2 + (m + m)^V 31 ) , (7.10) 

where /c and m are, respectively, the (1,1) and (1,2) terms of the Hcrmitian matrix Wl- 

Similarly, if we assume that V [0] is non-zero, then we can rearrange this formula to obtain [i and 
K in terms of the unknown v and the known tp\°\ tp^\ ip^ and ip^. We then obtain the following 
formulas for the trace and determinant of Wl — Wr: 

traced - Wr) = trace(^) + V ^ ^ J , (7.11) 

det(W L -W R ) = D + -^—(u^ - u^Wlu^), (7.12) 



where the constant D in (7.12) is given by 

D = det(Wi) + (n(^ [1 V [2] - - (^ [2] ) 2 + (m + m)^°ty p] ) • (7.13) 

Here n is the (2,2) term of Wl- 

We shall now prove a number of results which indicate how the Lj can be constructed to obtain 
spectral exactness. The following lemma shows that for a given A* below the essential spectrum, 
we can find a nearby value A so that det Uz,(bj, A) / and a(bj, A) = 1. Recall that hp h as been 
defined as the Friedrichs extension of the operator Lq, which has boundary conditions (7.1). 
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Lemma 7.1 Suppose that 

(a) The operator L is bounded below; 

(b) The interval [A*, A* + e] lies below the essential spectrum, and contains no eigenvalue of the 

Friedrichs extension Lp. 

Then 

(1) There exists /3q £ [a, b) such that for any (3 £ [/3q, b) and A £ [A*, A* + e], det Ul(/3, A) ^ 0; 

(2) If j3 £ [/?o, b) and u^,{(3) =fi 0, then there exists A £ [A*, A* +e] and a 2x2 real, symmetric matrix 

Wr such that: 

(i) v^((3) = W R u^(p); 

(ii) o~(/3, A) = 1 /or t/ie truncated eigenvalue problem 

£(y) = Ay, a < a < (3. 
Aiu(a) + A 2 v(a) = 0, 
v(l3) = Wr«(/3). 



Proof: (1) For a < (3 < b, consider the eigenvalue problem 

t(y) = ^y, a<x < (3, 

Aiu(a) + A 2 v(a) = 0, 
u{(3) = 0. 

We shall denote this eigenvalue problem by EP(/3). Let \ik denote the Friedrichs eigenvalues, and 



Hk{@) the eigenvalues of EP(/3). By Theorem 4.1 ^ fc (/3) \ as (3 / ' b 



Since A» + e lies below the essential spectrum, there are only finitely many Friedrichs eigenvalues 
(perhaps none) which are less than A* + e. 

If there are no Friedrichs eigenvalues below A* + e, then /ife(/3) > A* + e for all (3 € [a, 6) and all 
indices fc. In this case, we take /?o = a. 

If there are some Friedrichs eigenvalues below A* + e, let hn be the largest one. Since [A*, A* + e] 
contains no Friedrichs eigenvalues, /Ujv < A*. This implies that there is (3q £ [a, 6) such that 
Hn((3) < A* for all (3 £ [0o,b). Furthermore, if /i/v+i exists (below the essential spectrum) then 
//jv+i(/3) > A*jv+i > A* + e for all /3 £ (a, 6). If //jv+i does not exist below the essential spectrum, 
then /lijv+i(/3) does not lie below the essential spectrum for any (3 £ (a, 6). Thus no eigenvalue of 
EP(/3) lies in [A*, A*+e] for /3 £ [/9 , 6). But A is an eigenvalue of EP(/3) if and only if det U L {(3, A) = 0. 
This shows that det U L (J3, A) ^ for all A £ [A», A* + e] and £ [/3 , 6). 

(2) For/?£ [/3 ,6), and A£ [A*, A* +e], det C7 L (/3, A) ^ 0. Therefore W L (/3, A) = F L (/3, X)U L (f3, A)" 1 
is defined. For the eigenvalue problem in (2) (ii) , a({3,\) is the number of negative eigenvalues of 
W L (J3, A) - Wr. Now cx(f3, A) = 1 if 

det(W L {(3,\)-W R )<0. (7.14) 

The matrix Wr is constructed as indicated in the calculations preceding this lemma. (The construc- 
tion guarantees that v^(f3) — Wru^{(3).) Sin ce u^ {(3) ^ 0, det (Wl((3, A) — Wr) can be calculated 
by one of the formulas (|j) if {(3) + 0, or ( [712| ) if V [0] (/?) ^ 0. We can then force the inequality 



(7.14) by an appropriate choice of n or v, provided that 



.T 



^-u^W L {p,\)u^^Q, (7.15) 
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where u^p and v$ are evaluated at [3. 
We claim that the inequality ( 7.15| ) is satisfied for some A G [A*, A* + e]. This follows from the 



fact (shown in Greenberg jl(J) that Wp(f3, A) is a strictly decreasing matrix function of A in any 
A-interval containing no zeros of det Up{(3, A). Thus if u^v^ — u^Wl(P, A*)u^, = 0, and if A* is 



increased slighty to A»», then the inequality (7.15) will be satisfied for A = A*». (Here we have again 



used the fact that u^(fi) ^ 0.) □ 

Remarks. Suppose that L is a self-adjoint extension of Lq with domain D{L) = D(Lo) © Span (ip). 

(1) The function if> cannot have compact support. Thus the assumption u^([3) ^ (in part (2)(ii) 
of the preceding lemma) is satisfied on some sequence bj f b. 

(2) Suppose that there is a sequence bj / b such that u^(bj) = 0, for all j. Then the eigenvalues of 
L below the essential spectrum coincide with those of the Friedrichs extension Lp. These may be 
approximated by regular truncated problems with Dirichlct boundary conditions at the points bj. 

Theorem 7.1 Suppose that the operator L is bounded below, and has eigenvalues Ao,...,A n be- 
low the essential spectrum. Suppose also that the Friedrichs extension Lp of Lq has eigenvalues 
/io, • . . , /i n below the essential spectrum. [Note: we are not assuming that these are the only eigen- 
values below the essential spectrum.] Then the following inequalities hold between the eigenvalues of 
L and Lp: 

Ao<Mo, Aifc-i < A fc < k = l,...,n. (7-16) 

Proof: The inequality Xk < fj,k is well-known (see, e.g., Dunford and Schwartz [0, Problem D2, 
p. 1544]). For the inequality /J,k-i < A& we proceed by contradiction. Suppose that A& < and 



fix A* and e so that [A*, A* + e] C (Afe, /Ufe_i). By Lemma 7.1 we can choose a sequence of points bj 
converging to b and a set of associated boundary conditions at bj which give a (bj,^) = 1 for some 
G [A*, A* + e]. Let Hi(bj) be the eigenvalues of the truncated problem on [a,bj] with Dirichlet 



boundary conditions at x = bj. By Lemma 4.3 we know that for all sufficiently large j we have 



A* + e < }ik— i < A*fe-i(6j). Since A) is the nu mbe r of eigenvalues (J-i(bj) l ess than A, this 



implies that for A G [A*, A* + e], S^ibj^X) < k — 1 in (7.6). Since a(bj,v^') = 1, this means that 
N(bj, v^) — 5p{bj, f^') + a(bj, p^') < k. Thus the truncated operator Lj has at most k eigenvalues 
less than and therefore at most k eigenvalues less than A*. But the singular operat or L has 
k + 1 eigenvalues less than A*. This violates the spectral inclusion guaranteed by Theorem 3.1. □ 



If A„ < (j, n , then we can prove the following sharper version of Lemma 7.1 (2). 
Lemma 7.2 Suppose that 

(a) L is bounded below; 

(b) A„ < A* < fj, n . 

Then there exists 70 G [a, b) such that if j3 G ["fo,b) and u$(f3) ^ 0, there is a 2x2 real, symmetric 
matrix Wr such that 

(1) v^) = W R u^{f3); 

(2) er(/3, A*) = 1 for the truncated eigenvalue problem 

L{v) = ty, a <x < [3, 
Aiu{a) + A 2 v(a) = 0, 
v(f3) = W R u{[3). 
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Proof: By Lemma 7.1, there is [3q such that if (3 E [Po,b) then det Ul((3, A*) ^ 0. Therefore, for 
the eigenvalue problem in (2), cx(/3, A*) is the number of negative eigenvalues of Wl{(3, A*) — Wr, 
so <r((3, A*) = 1 if 

det(W L (P,\*)-W R )<0. (7.17) 



We construct Wr by the calculations preceding Lemma 7.1. The trace and determinant of 
W L (0,\*) ~ Wr are given by the formulas (^Tl}7_(|7^) , < ^UA^ , ( f7~l~2| ). The inequality ([TUT]) will 
be forced by an appropriate choice of re or v in (7.5) if 



(7.18) 



If the left hand side of (|7lg|) is zero, then det (Wi(/3, A*) - Wr) equals C = C(/3) in (fF7To| ) if 
V'' 1 '(/3) 7^ 0) or it equals Z? = £>(/?) in ( 7.13 ) if tp^(P) 0. If these coefficients are negative, then 
a((3, A*) = 1. To prove the existence of 70 6 [/?o, &) with the properties stated in the lemma, we shall 
argue by contradiction. Suppose that there is a sequence bj f b with the following properties: 

(i) u^(bj) ^ 0; 

(ii) u^bjfv^bj) - u-^(bj) T WL(bj, X)u^(bj) = 0; 

(iii) C(bj) > if V [1] (M ^ 0; D{b 3 ) > if W° ] (^) ^ 0. 

Note that C(bj) and D(bj) do not involve re or ^. Therefore det (W^(6j, A») — VFr) is independent 
of re or z^, and is nonnegative. By an appropriate choice of re or v, we can construct Wr so that 
trace {Whibj, A*) — Wi?) > 0. This implies that a(bj, A*) = for the eigenvalue problem 

t(y) = Ay, a<x<&,-, 
Aiu(a) + A 2 v(a) = 0, 
= W R u(bj). 

We shall denote this truncated eigenvalue problem by EP(&j). 

Let fik(bj) denote the eigenvalues of the truncated problem on [a,bj] with Dirichlet boundary con- 
ditions at x — bj. By Lemma 4.3, for sufficiently large j, 

A„ < A* < fi n < fi n (bj)- 

Since 5z,(6j,A*) is the number of eigenvalues fik(bj) less than A*, this implies that di,(bj,\*) < n, 
and 

N(bj, A*) = 5 L (bj,\*) +a(b j ,X») < n. 

Thus the approximating problems EP(6j) have at most n eigenvalues less than A», while the singular 
problem has n + 1 such eigenvalues. This contradicts spectral inclusion. □ 



Theorem 7.2 Suppose that the hypotheses of Theorem 1A_ hold, and suppose also that \ n < fi n . 
Then we can construct a sequence of operators Lj on intervals [a,bj] whose eigenvalues Xk(bj) 
converge, for each < k < n, to the corresponding eigenvalues Xk of L: 

lim Xk{bj) = Xk, k = 0, l,...,n. 

j— >oo 

Proof: Fix A* 6 (X n ,fi n ), and let bj /* b. By the previous lemma, we can construct truncated 
eigenvalue problems EP(&,): 

^(y) = Ay, a<x<bj, 
Aiu(a) + A 2 v(a) = 0, 
v(bi) = W jU (bj), 
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such that cr(6j,A*) = 1. Let Hk(bj) denote the eigenvalues of the truncated problem on [a,bj] with 
Dirichlet boundary conditions at x = bj. For sufficiently large j, 

Mn-l < f*n-l(bj) < K < Mn < Mn(^')- 

This implies that SlQjj, A*) = n, and 

N(bj,\*) =5 L (b ji )^)+a(b j ,K) =n + l. (7.19) 

Thus the approximating problems EP(&y) have n + 1 eigenvalues less than A*, and the same is true 
for the singular problem. Note that X n +i(bj) cannot converge to any of the singular eigenvalues less 



than A*, since equation (7.19) implies that A* < \ n +i(bj). The result now follows from spectral 



inclusion. □ 

This theorem shows that it is possible to construct a sequence of regular fourth order Sturm- 
Liouville problems to approximate the eigenvalues below the essential spectrum of a lim-3 singular 
problem, with each such eigenvalue being approximated by the corresponding regular eigenvalues 
of the same index. In other words, we can avoid having a situation in which, say, the eigenvalues 
Ao(6j) — > — oo as j — > oo. This would be an undesirable phenomenon since it essentially means 
that the regular approximating problems possess a spurious eigenvalue which is not approximating 
anything in the spectrum of the problem which interests us. However, the implementation of the 
theorem by a numerical procedure requires that we find a point A» in (A n , /i n ) for some n. We may 
not know enough about the spectrum in advance to be able to choose such a point. The following 
proposition can be useful for this purpose. 



Proposition 7.3 Suppose that the hypotheses of Theorem 7.1 hold. In particular, assume that 
A n < M«- 

(i) Suppose that for some k < n, Xk < fik- Then for each A £ (Xk,^k), the minimum value of 

o~(bj, A) is 1 for all sufficiently large j . 

(ii) Suppose that for some k <n, fJ-k-i < A&. Then for each A G (A*fc-i;Afc) the minimum value of 

o~(bj,\) is for all sufficiently large j. 

(iii) Suppose that A < Ao- Then the minimum value of o~{bj,X) is for all sufficiently large j. 
Proof: The proofs of the three parts of this theorem are all quite similar. Let Lj indicate the 



approximating truncated operators. From Corollary 3.1 we know that we can choose the Lj to 
achieve spectrally inclusive eigenvalue convergence. In Case (i), if A £ (A&, (ik) then Lj has at least 
fc+1 eigenvalues less than A for all sufficiently large j. The number of eigenvalues of Lj which are less 
than A is also given by ( |7.6| ) in which 5i,(bj, A) is exactly equal to k. This means that a~(bj, A) > 1. 
However we also showed above that we can choose the matrix Wr to achieve o~(bj, A) = 1, so this is 



the minimum value. For Case (ii), we exploit Theorem 7.2 to assert that we can choose the sequence 
Lj to be spectrally exact. In this case we know that for all sufficiently large j, each Lj will have 
precisely k eigenvalues less than A because \k-i < Mfc-i < A < A&. Again we know that for all 



sufficiently large j, we have 6i,(bj, A) = k for all sufficiently large j. Thus from (7.6), this spectrally 
exact sequence Lj must be giving us o~(bj,\) = for all sufficiently large j, which is clearly the 
minimum value that a can have. Finally, the proof of Case (iii) is virtually identical to the proof of 
Case (ii). □ 

This Lemma gives us a prescription for obtaining eigenvalue approximations. Suppose we want 
to compute, for each A below the essential spectrum of L, the number of eigenvalues of L which are 
less than A. Then all we need to do is set up a regular approximation over a truncated interval [a, bj], 



with bj sufficiently close to b, and choose our boundary condition (7.4) to minimize o~(bj, A) for this 



value of A. N(bj, A), which we can compute if we can solve a regular fourth-order Sturm-Liouvillc 
problem, will be the eigenvalue count which we seek. 
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8 Two singular endpoints 



8.1 The lim-2, lim-2 case 

In this case, for a problem posed over (a, 6), the preminimal domain is the set of functions with 
compact support in (a, b) and the minimal domain is its closure in the graph norm. However, in 
this case the minimal operator is self-adjoint. Thus, if it is bounded below, it is its own Friedrichs 
extension, and the results of Section [I] apply. If it is not bounded below then it is easy to see, by 



analogy with Proposition 6.2, that if there are gaps in the spectrum then a sequence of approximating 
truncated problems may be constructed in such a way that all will have an eigenvalue at some fixed 
point in a gap. Further information about the spectrum in such cases would therefore require 
approximations to the spectral function. 

8.2 The lim-2, lim-3 case 

We consider now a problem posed over an interval (a, b) in which x — a is a lim-2 endpoint and 
x = b is lim-3. At x — a we have no boundary conditions, while at x — b we have a single boundary 
condition of the form [y,ip](b) — 0. We may suppose that ip is real. We shall assume that the 
minimal operator is bounded below. 

If the boundary condition function ip yields the Friedrichs extension Lp then we have no more 
work to do: we have already seen how to obtain spectral exactness in this case. Other cases require 
a much more delicate treatment. We shall continue to denote the eigenvalues of Lp by fj,k] the 
eigenvalues subject to [y,ip](b) = will be dnoted by Afc. 

Theorem 8.1 Suppose that there is some n such that \i n lies strictly below any essential spectrum 
and suppose that X n < fi n . Let aj be a sequence such that aj \ a as j /* oo, and suppose that 
Dirichlet boundary conditions are imposed at x — aj. Then there exists a sequence bj with bj f b, 
and a set of boundary conditions at each bj including the boundary condition [y,i/>](bj) — 0, such 
that the eigenvalues A^ of the associated truncated eigenproblems satisfy 

lim a1 j) = A fe , k = 0,...,n. (8.1) 

Proof: Let Afc(a.,-, b) denote the fcth eigenvalue of the problem with one singular endpoint at x = b, 
with boundary conditions y(cLj) = = y'(aj), and [y,ip](b) = 0. Whatever point bj we choose and 
whatever boundary condition we impose there, we shall always have 

\X[ 3) - A fc | < |a£° - A fc ( ai ,6)| + \X k ( aj ,b) - A fe |. (8.2) 



Theorem 4.2 



gives 



lim |A fc (o 3 -,6)-Afc| = (8.3) 

j— >oo 

so it remains only to show that we can choose the bj and associated boundary conditions to get 

lim |A^ } -X k ( aj ,b)\ =0. (8.4) 

j->oo 

We do this by noting that we are now dealing with problems having just one singular lim-3 endpoint 
at x — b, the point x = aj being regarded as regu lar and fixed. Provided the appropriate hypotheses 



are satisfied, wc shall be able to apply Theorem 7.2 



We denote by fik (aj > b) the eigenvalues of the problems over (aj , b) with Dirich let conditions at 
x = q 7 a nd Friedrichs conditions at the singular endpoint x = b. From Theorem <L2 we have not 
only ( |8.3| ) but also 

lim \n k (a,j,b) - = 0. (8.5) 
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Thus the inequality A„ < \i n translates to X n (aj,b) < ii n (a,j,b) for all sufficiently large j. This allows 
us to apply Theorem 7.2 to deduce that there is a point bj < b and associated boundary conditions 



at x = bj, including the condition [y,if>](bj) = 0, such that 

\\®-\k( aj ,b)\ < ~. 



Combining this with (8.3) and (p.2|), our proof is complete. □ 



8.3 The lim-2, lim-4 case 

In this section we consider a problem posed over an interval (a, b) in which x — a is a lim-2 endpoint 
and x = b is lim-4. At x = a we have no boundary conditions. At x = b we shall impose two 
boundary conditions 

[y,u 3 }(b) = 0, [y,u 4 }(b)=0 (8.6) 

where, by recourse to the earlier results, we may assume that 113 and 1*4 are, for some real A, 
solutions of the differential equation which are square integrable at b, linearly independent relative 
to the minimal domain, and satisfy [1/3,143] = [1x3,144] = [1x4,144] = 0. As usual we denote the 
resulting self-adjoint operator by L and its eigenvalues by A^. 

We shall consider a sequence of approximating operators Lj defined over intervals (<Xj, bj) where 
Qj — ► a and bj — > b as j : 00. The boundary conditions (3.6) will now be replaced by 



y(aj) = = y'(a j ), [y, u 3 ](b 3 ) = = [y, u 4 ](b 3 ). (8.7) 

We shall denote the eigenvalues of these approximating problems by X k (aj,bj). 

Following the results which we presented in the case of just one lim-2 endpoint it is clear that 
for spectral exactness we had better not have spectrum extending to — 00 with a gap. With a lim-4 
endpoint we could, of course, have discrete spectrum extending to —00. We shall assume that this 
is not the case. Our main result is the following. 

Theorem 8.2 Suppose that L is bounded below and possesses at least n + 1 eigenvalues Ao, . . . , A n 
strictly below any essential spectrum. Suppose also that the sequence (aj)° c L 1 converges monotonically 
to a. Then we have the spectrally exact convergence 

lim Afc(aj, bj) = A&, fc = 0, . . . ,n. (8-8) 

j—*oo 

Proof: To establish this result we compare the following eigenvalues: 

• the eigenvalues Xk(aj,bj); 

• the eigenvalues \k(ai,bj) for problems on intervals (ai,6j); 

• the eigenvalues Afc(a, bj) of problems with a lim-2 singular end at x = a and a regular end at 
x = bj, with boundary conditions [y, u^](bj) = = [y, 144] {bj); 

• the eigenvalues Afc(a, b) := A& of L. 

We start with the triangle inequality, which gives 

|Afc(t4j,&j) - Afc(a, 6) < \\ k (a,j,bj) - \ k (a,bj)\ + \X k (a,bj) - X k (a, b)\ (8.9) 
Next, we observe that for all i < j, we have 

Xk{a,bj) < Xk{aj,bj) < \ k (a,i,bj). (8.10) 
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This is because a < aj < a; and the boundary conditions at the left hand endpoints are always 
Dirichlet. This inequality yields 

\Xk(dj,bj) - Afe(a, bj)\ < \Xk{ai,bj) - X k (a,bj)\ (8-11) 



Substituting ( p. 11 ) into the right hand side of ( |8.9| ) yields 

\\ k (a,j,bj) - X k (a,b)\ < |A fc (a, ; , bj) - Xk(a,bj)\ + \X k (a,bj) - X k (a,b)\ Vj > i. (8.12) 



If we now let j tend to infinity on the right hand side of ( 8.12 ) then the second term, \X k (a, bj) — 
X k (a, 6)|, will tend to zero by the results of Section ^, while the first term will tend to \X k (a,i, b) — 
Afc(a, b)\. Thus we have 

lira \X k (aj,bj) - Xk(a,b)\ <\X k (ai,b) - X k (a,b)\ Vi (8.13) 

j— >oo 



To complete the proof one need only show that the right hand side of ( 8.13 ) tends to zero as i — > oo. 



This is an immediate consequence of Theorem 4.2 



□ 



8.4 The lim-3, lim-3 case 

Since the case of one lim-3 endpoint was so awkward it should come as no surprise that the case of 
two lim-3 endpoints is the most difficult to treat. We shall consider an eigenvalue problem 

£y = Xy, xe(a,b), 

M](a) = 0, (8-14) 
[vM (b) = 0, 

which we denote by EP. We assume that the minimal operator associated with this problem is 
bounded below, that <f> and i\) are real, and that neither nor ip defines a Friedrichs boundary 
condition (since Friedrichs boundary conditions have been treated earlier). We shall denote the 
eigenvalues of EP by X k . We also require the problem 

£y = Xy, xe(a,b), 

[y,^}(a) = 0, (8.15) 
Friedrichs BC at b. 

which we denote by EP F and whose eigenvalues we denote by /x^, and the problem 

£y = Xy, xe(a,b), 

Friedrichs BC at a, (8.16) 
Friedrichs BC at b. 

which we denote by EP FF and whose eigenvalues we denote by [i k F ■ By variational methods it is 
easy to show that when the relevant eigenvalues all exist and lie below any essential spectrum, then 

Aft < fi k < fJ. k F - 

In line with our assumptions in Section 0, we shall assume that for some positive integer n, 

X n < fj, F < fi FF . (8.17) 



Our main result is that if, for some n, the inequalities in ( p. 17 ) hold, then we can obtain a sequence 
of approximating regular problems to give spectral exactness for the first n + 1 eigenvalues. 
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Theorem 8.3 Suppose that for some integer n, the eigenproblems EP, EP F and EP FF all have at 
least Ti + l eigenvalues strictly below the essential spectrum, and suppose that ( 8.1!\ ) holds. Then 
given any seguence (bj) with bj f b as j f oo, we can construct a sequence (aj) with aj \ a 
as j f oo, and a sequence of regular problems EP{a,j,bj) on the intervals [dj,bj] whose boundary 
conditions include [y, <j)](cij) = and [y,ijj](bj) = 0, such that for k = 0, . . . , n, the kth eigenvalue 
Xu of EP(aj , bj ) satisfies 

lim \® = A fc . 



Proof: Since ( S.17 ) holds, we can fix a point A* in (A„, /u, F ) such that A* fi F for all k. Now given 
the sequence bj consider the sequence of singular eigenproblems EP F (a, bj) defined by 



£y = Xy, xe(a,b), 

M](a)=0, 
y(bj) = - y'(bj), 

whose eigenvalues we denote H F (a, bj). From Section ^ we know that 

Mfe( a J 6 j)\Mfc asj/oo, k = 0,...,n. 

Similarly, we can set up the eigenproblems EP i?i? (a, bj) defined by 

£y = Xy, xe(a,b), 
Friedrichs BC at a, 
y(bj)=0 = y'(bj), 

whose eigenvalues we denote /i FF (a, bj). Again, Section || tells us that 

Mfc F (a, bj) \ ii FF as j / oo, k = 0, 



,n. 



(8.18) 



(8.19) 



(8.20) 



(8.21) 



Next, we examine EP F (a, bj) for each j. This is a problem with one lim-3 endpoint and one regular 
endpoint; we want to exploit the results of Section [t] to approximate its eigenvalues. To this end 
we must check the hypothesis of Theorem 7.2: do we have [i F (a,bj) < [i FF {a,bj)l Combining 
( [3.17 ), (3.19) and (3.21), it is clear that we do. Thus Theorem |7.2| can be applied: we can choose 
an endpoint at > a and a set of boundary conditions at on including the condition [y, 4>]{on) = 0, 
such that the eigenvalues, let us call them /j, F (ai,bj), of the resulting regular problem EP F (ai, bj) 
approximate the eigenvalues of EP F (a, bj): 



lim \fi F (ai,bj) - fi F (a, bj)\ =0, k = 0, . 



.22) 



Choose i — i(j) such that 
Set aj = ojj(j) . Then using 



\fi F [ai,bj) - fi F (a,bj)\ < ~. 



.IS) we have 



lim n F (a v bj) = fj, F , k = 0, 



.23) 



The problems EP F (aj, bj) will have the form 



ty = Ay, x e (a, b), 
v(aj) = W L (aj)u(aj), 
y(bj) = = y'(bj), 



(8.24) 
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for some matrices W^aj) chosen so that the boundary condition at x = a,j is satisfied by (j>. 

Given the problems EP (a,, bj) we now aim to change the boundary condition at x — bj. With 
Wi,(aj) fixed we attempt to choose a matrix Wn(bj) such that the eigenproblems EP(a,j,bj) defined 
by 

£y = Ay, x € (a,b), 

v( aj ) = W L { aj )u( aj ), (8.25) 
v{b 3 ) = W R (b 3 )u(bj), 

have precisely n + 1 eigenvalues less than A* . We must do this in such a way that the boundary 
condition at x = bj is satisfied by ip in order to guarantee spectral inclusion. Choose e > such 
that [A*, A* + e] C [A*,/j, f ). For each A £ [A*, A* + e] we can integrate the Hamiltonian form of 
the differential equation forward from x — a,j, starting with initial conditions C/j,(a,-,A) = I and 
Vt(aj,A) = Wl(oj), to obtain !7l(6j, A) and Vt(6j,A). 

We assert that if e is chosen sufficiently small then we shall have det Ui,{bj, A) ^ for all 
sufficiently large j and for all A £ [A* , A* + e] . For if this were not true then we could extract a 
subsequence of the bj tending to b and a sequence of values of A tending to A* at which we had 
det Ul = 0. These values of A would be eigenvalues of a sub seque nce of the problems EP F (aj,bj) 
and would therefore have to converge to /i F for some k, by (3.23). This would mean that we had 
A* = fi F for some k, which would contradict the choice of A *. 

If Utf,[bj) 7^ for sufficiently large j, then Lemma |7.2| implies that we can choose the matrix 
Wn(bj) so that the boundary condition v(bj) = Wuu(bj) is satisfied by if) and gives the eigenproblem 
EP(aj, bj) precisely one more eigenvalue below A* than the eigenproblem EP F (aj, bj): in other words, 
EP(cij, bj) has the same number of eigenvalues below A* as EP(a, b). Spectral inclusion (Lemma 3.4) 
now forces the convergence 

lim Air = Afc, k = 0, . ..,n. 

j— >oo 

If we have u^p(bj) = on a subsequence of the bj then, as in the case of one lim-3 endpoint, we can 
show that the eigenvalues \ J must converge to \x F for k = 0, . . . , n. By spectral inclusion, together 
with the ine quali ty A^ < fj, F , this then implies that A& = fi F for k — 0, . . . ,n, contradicting the 



assumption (8.17). This completes the proof. □ 



8.5 The lim-3, lim-4 case 

Suppose that x — a is a lim-3 endpoint, x — b is lim-4, and L m i n is bounded below. Let L be a 
self-adjoint extension of L m i n with boundary conditions [y, <f>](a) = 0, [y, ipi](b) = = [y, ip2](b). 
We may assume that <f> is real and ipi , ^2 are solutions of the differential equation with real A. We 
are interested in approximating the eigenvalues of L below the essential spectrum. Let Lq be the 
extension of L min with boundary conditions [y, ipi](b) = = [y, ^(b), and let M be the Friedrichs 
extension of L . Suppose that M has boundary conditions [y, 6](a) — 0, [y,^i](b) = = [y, il>2[(b). 
Let Ao < Ai < A2 < • • • be the eigenvalues of L below the essential spectrum, and let /iq < /^i < 
[i2 < • • • be those of M, Then A& < /ife. 

Theorem 8.4 Suppose that X n < /Lt n for some n. Then we can construct a sequence of regular 
operators Lj on truncated intervals [aj,bj], with eigenvalues such that 

lim Ai = A fe , k = 0, 1, . . . , n. 



Proof: Let a < bj < b, bj f b as j — > 00, and consider the following operators on the truncated 
interval (fl,bj\: Lj has boundary conditions [y,4>](a) — 0, [y, ipi](bj) = = [y, ip2}(bj)', Mj has 
boundary conditions [y, 9](a) — 0, [y,ipi](bj) = = \y,ip2[{bj). Lj is a truncation of L, and Mj is 
a truncation of M. By Corollary fH\ the sequences Lj and Mj are spectrally exact for L and M, 



33 



respectively. All of these operators have the same essential spectrum. Let X^ be the eigenvalues 
below the essential spectrum for Lj, and let /tL be those for Mj. Thus 

lim xi j) = A fc , lim 4 J) = Mfe- (8-26) 

j — * OO j — > OO 

Since A n < ^ n it follows that for sufficiently large j 

(8.27) 

Let L° be the operator on (a, &.,■] with boundary conditions [y,ipi](bj) = = [y, ip2](bj); then Mj is 



the Friedrichs extension of L°. Thus the hypothesis of Theorem 7^2 is satisfied for Lj. Consequently, 



for a sequence a, \ a, we can construct regular operators Ljj on intervals [cti, bj], with eigenvalues 



A,? '"^ , such that 



Km X ( *> j) = X { j ] , fc = 0,l,...,n. (8.28) 



For each j, choose i = such that 



l4 <J) -A^|<j, fc = 0,l,...,n. (8.29) 

Let aj = Qfi(j), £j = j, A^ 5 = A^* . Then is a regular operator on the interval [aj,6j] 
with eigenvalues XY ' , and lim XY = Afc, for fc — 0, 1, 2, . . . , n. □ 
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